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ABSTRACT 

Nearly all previous discussion on beaming effects in gamma-ray bursts have 
assumed a conical geometry. However, more and more observations on rela- 



> 

cn 

\l ' tivistic jets in radio galaxies, active galactic nuclei, and "microquasars" in the 

CN ; 

I Galaxy have shown that many of these outflows are not conical, but cylin- 

drical, i.e., they maintain constant cross sections at large scales. Thus it is 
■ necessary to discuss the possibility that gamma-ray bursts may be due to 

I ' highly coUimated cylindrical jets, not conical ones. Here we study the dynam- 

o ■ 

|H I ical evolution of cylindrical jets and discuss their afterglows. Both analytical 

, and numerical results are presented. It is shown that when the lateral ex- 

pansion is not taken into account, a cylindrical jet typically remains to be 
• ■ 

' highly relativistic for ^ 10^ — 10^ s. During this relativistic phase, the opti- 

$^ ' 

' cal afterglow decays as S^, oc t~P^'^ at first, where p is the index characterizing 

the power-law energy distribution of electrons. Then the light curve steepens 
to be Siy cx t~(P+^^/^ due to cooling of electrons. After entering the non- 
relativistic phase (i.e., t > 10" s), the afterglow is cx i-(5p-4)/6 
if the cylindrical jet expands laterally at co-moving sound speed, then the 
decay becomes cx t'P and cx f-(i5p-2i)/io _ ^-(i5p-20)/io ^^le ultra- 
relativistic and non-relativistic phase respectively. Note that in both cases, 
the light curve turns fiatter after the relativistic-Newtonian transition point, 
which differs markedly from the behaviour of a conical jet. It is suggested that 
some gamma-ray bursts with afterglows decaying as t~^'^ — t~^ '^ may be due 
to cylindrical jets, not necessarily isotropic fireballs. 
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1 INTRODUCTION 

The discovery of gamma-ray burst (GRB) afterglows in early 1997 (van Paradijs et al. 1997; 
Sahu et al. 1997) has revolutionized the researches in the field (Wijers 1998). Since then, 
X-ray afterglows (Costa et al. 1997; Piro et al. 1998) have been detected from more than 20 
GRBs, of which a half were observed optically (e.g., Groot et al. 1997; Galama et al. 1997; 
Metzgcr ct al. 1997; Kulkarni et al. 1998; Pedersen et al. 1998; Garcia ct al. 1998; Akerlof et 
al. 1999; Vreeswijk et al. 1999) and one third were even observed in radio bands (e.g., Frail 
et al. 1997; Galama et al. 1998a, b). The cosmological origin is thus firmly established, and 
the so called fireball model (Goodman 1986; Paczyhski 1986; Rees & Meszaros 1992, 1994; 
Meszaros & Rees 1992; Meszaros, Laguna & Rees 1993; Meszaros, Rees & Papathanassiou 
1994; Katz 1994; Sari, Narayan & Piran 1996) becomes the most popular model, which is 
found successful in explaining the global features of GRB afterglows (Meszaros & Rees 1997; 
Waxman 1997a, b; Tavani 1997; Vietri 1997; Wijers, Rees & Meszaros 1997; Meszaros, Rees 
& Wijers 1998; Sari, Piran & Narayan 1998; Dermer et al. 1999a, b; Wijers & Galama 1999; 
Dai & Lu 1998, 1999). For detailed reviews on recent progresses, see Piran (1999) and van 
Paradijs, Kouvehotou & Wijers (2000). 

However, we are still far from resolving the puzzle of GRBs, because many crucial prob- 
lems arc largely open. For example, whether a GRB is due to a highly collimated jet or an 
isotropic fireball is still controversial. This issue has been discussed extensively. It is generally 
believed that due to both the edge effect (Panaitescu & Meszaros 1999; Kulkarni et al. 1999; 
Meszaros & Rees 1999) and the lateral expansion effect (Rhoads 1997, 1999), afterglows 
from a jet will be characterized by a break in the light curve during the relativistic phase. 
GRBs 990123, 990510 and 000301c are regarded as good examples (Kulkarni et al. 1999; 
Harrison et al. 1999; Castro- Tirado et al. 1999; Sari, Piran & Halpern 1999; Wijers et al. 
1999; Rhoads & Fruchter 2001; Bcrgcr ct al. 2001; Jensen et al. 2001; Masetti et al. 2000). 
But detailed numerical studies show that the break is usually quite smooth (Panaitescu & 
Meszaros 1998; Moderski, Sikora & Bulik 2000). Recently Huang et al. (1998a, b, 1999a, 
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b) stressed the importance of the non-relativistic phase of GRB afterglows. They have used 
their refined dynamical model to discuss afterglows from jetted GRB remnants (Huang et al. 
2000a, b, c). It is found that the most obvious light curve break does not appear in the rel- 
ativistic phase, but occurs at the relativistic-Newtonian transition point. They also stressed 
that the break is parameter-dependent (Huang et al. 2000b). Based on their investigations, 
the rapid fading of optical afterglows from GRBs 970228, 980326, 980519, 990123, 990510 
and 991208 has been suggested as evidence for beaming in these cases (Huang et al. 2000d). 

In all these discussion, a conical jet (i.e., a jet with a constant half opening angle; of 
course, in some more realistic discussion, the conical jet is permitted to expand laterally) 
has been assumed. However, we should note that this might not be the case. For example, we 
know that a relativistic flow (i.e., a jet) is a common feature of radio galaxies. Interestingly, 
copious examples have shown that these outflows are not conical, instead, they maintain con- 
stant cross sections at large scales (Perley, Bridle & Willis 1984; Biretta, Sparks & Macchetto 
1999; also see GLAST home page: \http://www-glast. sonoma.edu/gru/agn/Default.htni and 
Chandra home page: ^ttp://chandra. harvard. edu/photo/cyclel/pictor/index.htni^ for good 
examples). Here we call a jet maintaining a constant cross section a cylindrical jet. It is 
reasonable to deduce that ultra-relativistic outflows in GRBs might also be cylindrical jets, 
not conical ones. In fact, this idea has already been suggested as GRB trigger mechanism 
by Dar et al. (Shaviv & Dar 1995; Dar 1998, 1999). 

Afterglows from a cylindrical jet is likely to differ from those of a conical jet markedly. 
If GRBs are really due to cylindrical jets, then we need to examine their dynamics and 
afterglows in detail, which is just the purpose of this article. The structure of our article is 
as follows. We describe our model in Section 2 and present the analytic solution in Section 
3. Section 4 is the numerical results and Section 5 is brief discussion. 



2 MODEL 

Let us consider a highly coUimated GRB remnant, ploughing its way through a homogeneous 
interstellar medium (ISM) where the number density is n. Denote the lateral radius of the 
remnant as a, and its distance from the central engine of the GRB as R. The remnant will 
maintain a constant cross section (i.e., a = const), except that it may expand laterally at a 
speed v_i_. 
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2.1 Dynamics 

In the literature, the following differential equation is usually used to depict the expansion 
of GRB remnants (e.g., Chiang & Dermer 1999; Piran 1999), 

^ = -^ (1) 
dm M ' ^ ^ 

where m is the rest mass of the swept-up medium, 7 is the bulk Lorentz factor and M 
is the total mass in the co-moving frame, including internal energy. However, Huang et al. 
(1999a, b) pointed out that in the adiabatic case equation (0) does not match with the Sedov 
solution during the non-relativistic phase. They have proposed a refined equation, 

±L = (2) 

dm Mej + em + 2(1 - e)7m' 

which is shown to be correct for both adiabatic and highly radiative shocks, and in both 
relativistic and Newtonian phase. Here Mej is the initial ejecta mass and e is the radiative 
efficiency. In this paper we adopt equation (^ to describe the dynamical behavior of our 
cylindrical jet. For simplicity, we will only consider the adiabatic case (i.e., e = 0). 

The evolution of radius (R) , the swept-up mass (m) , and the lateral radius (a) is described 
by (cf. Huang et al. 2000a, b), 

^ = /5c7(7 + v'7'-l), (3) 
dm 2 

— = na nmp, (4) 
f = ^±(7 + v/^), (5) 



where t is the observer's time, (3 = 1/7^ — I/7, c is the speed of light, and m^ is the mass of 
proton. As for v±, we have two possibilities, just like in the case of a conical jet: (i) v± = 0, 
this means the jet does not expand laterally, it maintains a constant cross section strictly. 
The problem is then greatly simplified; (ii) v± = Cg, i.e., the jet expands laterally at the 
co-moving sound speed, of which a reasonable expression has been derived by Huang et al. 
(2000a, b) as, 

= 7(7 - 1)(7 - ^)-r-4 TT^'' 

1 +7(7 - 1) 

where 7 ~ (47 + l)/(37) is the adiabatic index (Dai, Huang & Lu 1999). In the ultra- 
relativistic limit (7 ^ 1,7 ~ 4/3), equation (|]) gives = c^/3; and in the non-relativistic 
limit (7 ~ 1,7 ^ 5/3), we simply get = 5/?^c^/9. 
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Equations (0) — @ give a thorough description of the dynamical evolution of our 
cylindrical jets. We will solve them analytically in the next section and then present our 
numerical results in Section 4. 

2.2 Synchrotron radiation 

A strong blastwave will be produced due to the interaction of the remnant with the ISM. 
Synchrotron radiation from the shock-accelerated ISM electrons gives birth to GRB after- 
glows. As usual, we assume that the magnetic energy density in the co-moving frame is a 
fraction of the total thermal energy density (i^'^/Svr = ^b^')) that electrons carry a 
fraction of the proton energy, which means that the minimum Lorentz factor of the ran- 
dom motion of electrons in the co-moving frame is 7e,min = Ce(7 — l)mp/me + 1, where rrie is 
the electron mass. In the absence of radiation loss, the distribution of the shock-accelerated 
electrons behind the blastwave can be characterized by a power-law function of electron 
energy, dN^/d'je oc %^ (7e,min < 7e < 7e,max), where 7e,max is the maximum Lorentz factor, 
7e,max = 10® (5'/! G) "^/^ (Meszaros et al. 1993; Vietri 1997; Dai et al. 1999). 

However, radiation loss may play an important role in the process. Sari et al. (1998) have 
derived an equation for the critical electron Lorentz factor, 7c, above which synchrotron 
radiation is significant: 

7c = Gnniec/iaT'yB'h), (7) 

where is the Thompson cross section. The distribution of radiative electrons (i.e., those 
with 7e > 7c) becomes another power-law function of dN^/d'je oc 'J^^^'^^^- A good description 
of electron distribution allowing for radiation loss has been presented by Dai et al. (1999, also 
see: Huang et al. 2000a, b). We will use their expression in this research. For convenience, 
let us define a characteristic frequency of Uc = 'j'j'^eB' / {2nmec) here (Sari et al. 1998). 

In the co-moving frame, synchrotron radiation power at frequency u' from electrons is 
given by (Rybicki & Lightman 1979), 

^ vs^ ...... ^ 

where e is electron charge, z/^j, = / {Anmec) , and 

F{x) = X / K,/sik)dk, (9) 

Jx 

with K^i2,{k) being the Bessel function. We assume that this power is radiated isotropically. 
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— = ^ . (10) 

Let 6 be the angle between the axis of the jet and the hne of sight, and define /x = cos B, 
we can derive the angular distribution of power in the observer's frame (Rybicki & Lightman 
1979), 

dP{u) 1 dP'{u') 1 P'(z/') 



dn 73(l-/5/i)3 dn' 73(l-/5/i)3 47r 



V 



(11) 

(12) 



7(1 -/i/5)- 

Then the observed flux density at frequency v is (Huang et al. 2000a, b), 

4 (^^) - iH^^^J' <^<^ - '^^"^ • '''' 

where A is the area of our detector and is the luminosity distance. 



3 ANALYTIC SOLUTIONS 

Since the geometry concerned here is still quite simple, we can derive approximate solutions 
to the problem analytically. In this section, for simplicity, we assume G = 0, i.e., the jet is 
just moving along the line of sight. Obviously, we then do not need to consider the effect of 
the equal arrival time surfaces (Waxman 1997b; Sari 1997; Panaitescu & Meszaros 1998). 
We discuss the two cases of f ^ = and v± = Cg separately. Note that we take e = all 
through the article. 

3.1 t>^ = case 

During the ultra-relativistic phase, since 7 1, /3 ~ 1, and 7m ^ Mgj, equations (||) — (H) 
reduce to, 

d^ 7 dm dR ^ n ,^ ^. 

!^ — , — — ^ const, — ^ 27^c. (14) 

dm 2m dR dt 

Then it is very easy to get: 

7 0cm~^/^ mocR, Roct^^^. (15) 

Additionally, from the shock conditions (Blandford & McKee 1976; Huang et al. 1998a, b), 
the co-moving thermal energy density is e' oc 7^ oc so we have B' 

and z/c oc 77^-8' oc t~^. Since the peak frequency of synchrotron radiation is oc 77emin-^' 
(Sari et al. 1998), we have oc t~^, and the peak flux density is Su,ma.^ oc Ne'jB'/^^'^ oc 
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iVe7^5' oc The factor of 7 ^ in S'jy^max is due to the relativistic beaming effect. At last, 

we can derive the flux density of the afterglow as follows, 

' ^.,max(l^)"^'"'^^'0CZ/-(*'-l)/2t-f/2, (z/^ < Z/ < sloW COoliug), 



S,, 



S,,^^,{^y^'"'^^\^y''^' cxiy-P/H-(P+^y^ {v^<v,<v, slow cooling). 



(16) 



Sua 



-1/2 



OC V ^, (z/f. < Z/ < I/rn, fast COollUG 



^ ^.,max(^) '^'(^) OC U-P/H-^P+'y\ {l^c < l^m < 1^, fast COoliug). 

In the non-relativistic phase, 7~l,/3-Cl, (i7~ pdp, equations (H) — (||) reduce to. 



-p/2 



dp (3 dm dR ^ 

- — ~ , — — ~ const, — — ~ pc. 

dm 2m dR dt 

The scaling laws become, 
/3ocm-^/^ mocR, R(xt^'^. 



(17) 



(18) 



According to the shock conditions (Huang et al. 1998a, b), the co-moving thermal energy 
density is e' oc 7 — 1 oc t^"^^'^. Then we have B' oc (e')^/^ oc t"^/^, Uc oc 77c-B' oc t^^, 
oc 77emin-^' t^^^^, S^^max OC 'jN^^B' OC t^^^ , and finally the flux density is, 

' ^.,max(l^)"^''"'^^'(XZ/-(*^-l)/2H5f-7)/6, (u^^ < U < U,, sloW COoliug), 



Sua 



-(p-l)/2 / , 



-1/2 



(^) OC zy-p/2t-(5p-4)/6^ <V^<V, sloW COollug) 



OC u 

-1/2^-1/6 

(z/c < < fast cooling) 



(19) 



^ >^.,maxfe) OC Z/-P/2r(5p-4)/6^ {^c < < fast COoliug). 

Let us discuss the decay of optical afterglows (e.g., R band, at the frequency z/r) in more 
detail. In the ultra-relativistic phase, t'm ^ '^r ^ is typically satisfied initially, then we 
have S'r oc z/"^^"^)/^^"^/^ at early stages. Since z/m and Vc both scale as t~^, we will eventually 
have < < then the afterglow at later stages is just oc z/~P/^t~*^P+^)/^, i.e., the 
light curve steepens by t^/^ due to the cooling effect of electrons. After entering the non- 
relativistic phase, since Ura oc t~^/^ and Uc oc t~^, the relation of z/m < i^c < t^R will continue 
to be satisfied. The afterglow then decays as Sr oc i'^p/'^t^(^P~^y^ ^ which means the light 
curve usually becomes flatter again, by t('^~2p)/6 



-p/2 
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3.2 = Cs case 

During the ultra-relativistic phase, 7 ^ 1, /5 ~ 1, 7m ^ Mej, Cg ~ c/^/3, equations 
(|) reduce to, 

dj 7 (im 2 '^'^ 27c dR ^ 2 

(im 2m' (ii? (it 1/3' dt 

The corresponding solution is, 

7 oc a oc t^^'^, m oc t. 



(20) 



(21) 



Similarly we have e' oc 7^ oc t \ 5' oc (e')^/^ oc t z/c oc 77^-6' oc oc 77c,min-^' oc t ^, 

and S'j.^max ^cT-B'/t^^ oc t^^. The flux density of the afterglow is then given by. 



S,, 



Su,ms.^ (^) ^^^^ OC 1/ (l^m < < ^^c, sloW COohug) 



-(p-l)/2 / , 



-1/2 



OC ^/^t (z/m < t'c < slow cooling) 

5'j.,max (^) ^ OC Z/^^/^t^^ {l^c < < t'm, fast COoliug) 



(22) 



^ ^i.,max (7^) (^) CX Z/ P/^^-^, (z/c < Z/^ < Z/, fast COoliug). 

Note that in most cases, the flux density decays as 

In the non-relativistic phase, 7 ~ 1, /3 <^ 1, Cg ~ v^/5c/3, equations (|) — (H) reduce 

to, 

P dm 2 

- — pa , — — vra nrrtp, — ~ 

dm 2m dR dt 3 

The solution is also easy to get as, 
/3 oc m"^/^ m(xR^, Rocaoct^^^. 



da y/bdc dR 



(23) 



(24) 

Then we have: e' oc 7 - 1 oc t~^/^, B' oc (e')^/^ oc z/c oc 77^-8' oc t"^/^, oc 77e,min-S' oc 

oc ^^,max OC 7A^e5' OC t^/s. So the flux density is, 

' 5.,n.ax(^)'^'"'^^'oCZ.-(^-l)/2t-(15^'-21)/10, (z.^ < < Z/„ sloW COoliug), 



-b l)/2 p/2 ^ ^_p/2^_(i5p-20)/10^ < Z/, < Z/, sloW COoliug) 



5.^ 



(25) 



Sua 



-1/2 



OC z/ 

-1/2^1/2 

(z/c < z/ < z^in, fast cooling). 



5.,max(^) '^'(^) "'^ ^u-P/H-^^^P-^'')n\ {v,<Vr.<y. fast COoliug) 
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Again let us have a look at the R band optical light curve. In the ultra-relativistic 
phase, typically we have <^ z^r -C Uc, so the afterglow decays as S'r oc i/~^P~^^^'^t~P. 
After entering the non-relativistic phase, we will have t'm < < t'c at first, which means 
S'r oc ^y-(p-l)/2i^-(l5p-2l)/lo^ ^Yiat the light curve becomes flatter by if:(2i-5p)/io_ gj^j.^ 

now decreases as Uc oc t^^^^, at sufficiently late stages, we will eventually have < Uc < z/r, 
which leads to S'r oc z/"P/2t-(i5p-20)/io^ j -(.j^g light curve is still flatter than that in the 
relativistic phase (by t(20-5p)/i0). in short, after entering the non-relativistic phase, the light 
curve is expected to become flatter obviously, which is quite different from the behaviour of 
a laterally expanding conical jet (Huang et al. 2000a, b, c, d). 

4 NUMERICAL RESULTS 

Now we solve equations (0) — (§) numerically and compare the optical (R band) afterglow 
light curves with those of conical jets. First we describe our initial values and parameters 
briefly. The ISM density (n), the original ejecta mass (Mgj) and the initial value of 7 are 
fixedly set as n = 1 cm^^, Mgj = 5 x IO^^^Mq, and 70 = 300, respectively. We further 
assume = 1.2 x 10^^ cm, then the initial values of m and R are given by: mo = Mgj/SOO, 
Rq = itlq/ {-Ko^nm^) ^ 1.2 x lO^'' cm. These values indicate that the angular radius of our 
GRB ejecta from the central engine is only ~ 10~^ at the beginning, and in the ultra- 
relativistic phase considered here it is just as powerful as an isotropic fireball with kinetic 
energy -Eiso,o ~ 10^"^ ergs. In calculating the R band fiux densities, we take Dl = 1 Gpc, and 
will assume p = 2.5, = 0.1, = 10~^, 9 = unless declared explicitly. 

Here, we would like to point out that many of our initial values and parameters (such 
as n, Mej, 7o, ao and Dl) have been arbitrarily chosen. However, variation of them within 
reasonable ranges will not alter the major conclusions we draw below in this section. 

4.1 v± = case 

Assuming no lateral expansion, we have evaluated equations (0) — (^ numerically. Fig. 1 
shows the evolution of the Lorentz factor, 7. The jet becomes non-relativistic after ~ 10^^ s. 
In the relativistic phase the slope is approximately —0.26 and in the non-relativistic phase it 
is about —0.68, both consistent with the theoretical value (i.e., —1/4 and —2/3 respectively, 
see Section 3.1). For comparison, the dashed line shows the evolution of 7 for a conical jet 
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(parameters: -Ejso.o = 10^^ ergs, 70 = 300, with an initial half opening angle 6*0 = 0.1). We 
see that the cylindrical jet decelerates much slower, which is easy to understand. 

Fig. 2 illustrates the Lorentz factor as a function of R. The slope of the solid hne is 
—1.00 in the relativistic phase and —0.50 in the non-relativistic phase, also consistent with 
the theoretical results. This further confirms the correctness of our programs. 

The R band afterglow light curves are shown in Fig. 3. The sohd line represents a 
cylindrical jet with p = 2.5. Its slope is approximately —1.26 before 10^ s, consistent with 
the theoretical value of —1.25 (for i^m < i^r < Vc case, see Section 3.1). The slope is —1.89 
between 10^° s and 10^^ s, also in good agreement with our analytical results (—1.75 for 
i^m < i^c < ^^r)- After ~ 10^^ s, the jet enters the non-relativistic phase and the light 
curve becomes flatter shghtly, just as required by the theoretical result. The dotted hne 
here represents afterglows from a conical jet (p = 2.5). We can see that the light curve of 
a cylindrical jet differs markedly from that of a conical one, mainly in the following two 
aspects: (i) the relativistic phase lasts for a very long period (> 10^'^ s typically), during 
which only a slight break (by i^/^, due to cooling of electrons) can be observed; (ii) after 
entering the Newtonian phase, it becomes flatter shghtly, not steeper. Please note that at 
the end point of each line, the characteristic Lorentz factor of electrons has already been as 
small as 7e ~ 5. 

Fig. 4 illustrates the effect of three parameters, ^'^'^ on the light curve. It is 

clear that does not change the shape of the light curve markedly. However, because z/c 
strongly depends on B' (7c oc v^. oc B'~^), we see that affects the cooling break 

seriously. In our calculations, since we take 70 = 300, we have found that for any © smaller 
than 1/300, the light curve will be almost identical to the © = case. But in the © = 0.1 
case, the flux density is lowered down markedly before ~ 10^ s, which means a cylindrical 
jet is observable only within a narrow solid angle. 

4.2 v± = Cg CEise 

When lateral expansion is included, the results are quite different. The solid line in Fig. 5 
illustrates the evolution of 7 in this case. We can see that the jet becomes Newtonian at 

~ 10^ — 10^ s, much earlier than the v± = case. Note that its slope is approximately 
—0.51 in the ultra-relativistic phase and —1.20 in the non-relativistic phase, both consistent 
with the theoretical values (i.e., —1/2 and —6/5 respectively, see Section 3.2). Fig. 5 also 



© 2001 RAS, MNRAS 000, 000-000 



11 

shows that when = Cs, a cyhndrical jet decelerates faster than a conical one, which is 
just contrary to the v± = case. The 7 - i? relation is shown in Fig. 6 by the solid line. 
In the non-relativistic phase, the slope is ~ —3.0, consistent with the theoretical value. 
However, in the ultra-relativistic phase, the 7 - it! relation is not a power-law function, but 
an exponential one. 

Fig. 7 illustrates the evolution of the jet profile on the y - z plane. Generally speaking, 
the lateral expansion is not significant during the ultra-relativistic phase. But in the non- 
relativistic phase, the lateral radius increases synchronically with the radial scale. 

R band afterglows from cylindrical jets expanding laterally are shown in Fig. 8. The solid 
line is plot with p — 2.5. Its slope is ~ —2.63 before 10^ s, consistent with the theoretical 
value of —2.5 (see Section 3.2). It turns flatter in the non-relativistic phase, also consistent 
with our analytical results. The dotted line corresponds to a laterally expanding conical 
jet {p = 2.5), where a break due to the relativistic- Newtonian transition can be clearly 
seen (Huang et al. 2000a, b, c, d). It is clear that the light curve of a cylindrical jet differs 
from that of a conical one mainly in the following aspects: (i) no obvious breaks could be 
observed; (ii) in the relativistic phase, the decay is always as quick as t'^, but for a conical 
jet the decay is only early stages (Rhoads 1997, 1999); (iii) after entering the 

non-relativistic phase, the light curve becomes flatter shghtly, also quite different from that 
of a conical jet. Fig. 9 illustrates the effect of .^c, © on the light curve. Similar to 

the v± = case, the flux density is markedly lowered down by a non-zero G, but generally 
speaking, and ^| do not change the slope of the light curve here. 

5 DISCUSSION AND CONCLUSIONS 

The importance of coUimation in GRBs goes beyond influencing the detailed shape of the 
light curve: it reveals something about the central engine, and affects our estimates of how 
many GRB progenitors there need to be (van Paradijs et al. 2000; Wijers et al. 1998; Lee, 
Brown & Wijers 2000). Most theorists now expect GRBs to be signiflcantly coUimated, 
because jets seem to be produced by all systems where matter is undergoing disk accretion 
onto a compact central object, including young protostars, "microquasars" (i.e., black hole 
binaries in the Galaxy), capture of a main sequence star by a moderate massive black hole in 
normal galaxies (Cheng & Lu 2000), radio galaxies, and active galactic nuclei (Lamb 2001). 
Observations on these systems have interestingly indicated that an astrophysical jet usually 
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maintains a constant cross section at large scales (Perley et al. 1984; Biretta et al. 1999). In 
fact, it has already been pointed out by Dar et al. (Shaviv & Dar 1995; Dar 1998, 1999) that 
although a conical beam can solve the "energy crisis" of GRBs, it suffers from some other 
deficiencies of isotropic fireballs, while a cylindrical jet does not. However, nearly all previous 
discussion on afterglows from coUimated GRB remnants have assumed a conical geometry, 
which, we believe is not enough. Here we have discussed afterglows from cylindrical jets in 
detail. Both analytic and numerical results arc presented. 

For a cylindrical jet without lateral expansion, in the relativistic phase the flux density 
typically decays as S^, oc i"^/^ at first. Then the hght curve steepens by t^/'^ due to cooling 
of electrons. After entering the Newtonian phase, the dechne becomes shghtly slower. Si, oc 
^-(5p-4)/6^ just contrary to the behavior of a conical jet. Note that the remnant usually 
decelerates very slowly in this case, i.e., 7 oc t~^/^ when 7 ^ 1, so that it keeps to be highly 
relativistic as long as t ~ 10* — 10^ s. When lateral expansion is included and assuming 

= Cg, the afterglow light curve is S^, oc in the ultra-relativistic phase, and it flattens 
to be S,y oc i-(i5p-2i)/io 0^ ^-(i5p-2o)/io ^Yie the non-relativistic phase, also quite 

different to a conical jet. 

In previous sections, we have assumed that the lateral expansion velocity of our cylindri- 
cal jet can be either v^ = Qorv^ = Cg and we have not given a preference for either of these 
two cases. Now we discuss this question. At the first glance, the v± = Cg seems more realistic. 
However, observations on relativistic outfiows in radio galaxies reveal that lateral expansion 
is usually unobtrusive (Perley et al. 1984; Biretta et al. 1999). For example, observed profiles 
are generally quite different from the profile in our Fig. 7. The high degree of coUimation 
and the suppression of lateral expansion is probably due to magnetic confinement (Dar 1998; 
Mirabel & Rodriguez 1999) and/or the external pressure. We believe that in many cases, 
even if the cylindrical jet expands laterally, its velocity will be v± ^ Cg. 

Afterglows from some GRBs, such as GRBs 970508, 971214, 980329 and 980703 have 
been observed to decay slowly and steadily within a long period (i.e., up to several months, 
see van Paradijs et al. 2000), with S^, oc t~^-^ — t~^-^. Although an isotropic fireball is 
the most natural and popular model for these GRBs, Huang et al. (1998a; also see Dai et 
al. 1999) pointed out that there is a problem: usually a typical fireball with kinetic energy 
-£"180,0 = 10^^ ergs will become non-relativistic in one or two months, then the theoretical light 
curve will become steeper slightly (Wijers et al. 1997; Huang et al. 1998b, 1999a, b). So it is 
difficult to understand why the light curve is a single straight line that lasts for 0.5 — 1 year 
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(e.g., GRB 970508). Here we suggest that this kind of GRBs might be due to cyhndrical jets 

(with v± <^ Cs): a hght curve with slope —1.1 1.3 can be easily produced by assuming 

2.2 < p < 2.6, and the steady dechne is explained by the fact that these cyhndrical jets 
are highly relativistic even when t ~ 10^ — 10^ s. Note that in this model, 7-rays of the 
initial burst are beamed into a solid angle with the angular radius ~ I/to- For To ~ 300, the 
beaming angle is ~ 0.003; and for 70 ~ 100, it is ~ 0.01, which means the energy involved 
will be greatly reduced. 

GRBs 990123, 990510, 991216 and 000301c are characterized by an obvious break in the 
afterglow light curve. A break steeper than i^/^ could not be explained directly by our simple 
cylindrical jet model, either v± = or v± = Cg. However, there are many other factors that 
have not been taken into account in this article, such as the inhomogeneity of ISM (especially 
n{R) oc R""^, Dai & Lu 1998; Chevalier & Li 2000), and the nature of the lateral expansion. 
For example, we can imagine that in some cases, since the shock is highly radiative at first 
(which means the remnant is relatively cool), the magnetic field is much stronger initially, 
and at the same time the external pressure may be high, then the lateral expansion can 
be successfully suppressed, which will lead to a flat light curve of Si, oc t"^ — t~^-^ at 
early stages. At later stages the lateral expansion takes effect and the light curve becomes 
Si, oc t"'^ — t~^. So an obvious break appears. In short, we suggest that the cyhndrical jet 
model should be paid attention to and should be investigated further since it seems to be 
more realistic than the conical jet model on the observational basis. 

Recently, Wang & Loeb (2001) discussed the emission from bow shocks of beamed GRBs. 
They found that the emitted flux from this bow shock may dominate over the direct emission 
from the conical jet for lines of sight which are outside the angular radius of the jet emission. 
The event rate for these hues of sight is thus greatly increased. For typical GRB parameters, 
they found that the bow shock emission from a jet of half-angle ~ 5° is visible out to tens 
of Mpc in the radio and hundreds of Mpc in the X-rays. In most of our calculations where 
6 = is assumed, the bow shock emission generally will not play a dominant role. However, 
we should bear in mind that the bow shock emission is obviously more important for a 
cylindrical jet than for a conventional conical jet, since observers of a cylindrical jet are 
more likely to be situated outside the initial emission cone. For example the bow shock 
emission may take effect in those © = 0.1 cases in Figs. 4 and 9. A direct effect of bow 
shock emission from cylindrical jets is that the beaming factor will be considerably reduced, 
so that it will be much easier for us to observe them. We will not go further to calculate 
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the bow shock emission here. But in the future, if the so called "burstless afterglows" are 
observed, then the question of bow shock emission will become very important. 

Finally we would like to say a few words about the origin of the cylindrical jet. It should 
not be born a cylinder. In fact, the outflow in a black hole-accreting disk system is likely 
to be conical initially, with the half opening angle ^ <^ 1. Assuming 9 ~ 0.1, then after 
propagating to a distance of i? ~ 10^^ cm, it will acquire a lateral radius of a ~ 10^^ cm. 
So in our calculation we have implicitly assumed that the jet has become cylindrical before 
R ~ 10^^ — 10^^ cm. We will not go to any further details of the mechanism that makes a 
conical jet into a cylindrical one here, they are beyond the purpose of this article. However, 
if our cylindrical jet is really originated from a conical outflow very near to the central black 
hole, then we can draw an interesting conclusion as follows. Adopting the parameters used in 
Section 4, our cylindrical jet has a total kinetic energy of E^j = 7oMejC^ 3 x 10^^ ergs. Since 
this energy comes from a cone with ~ 0.1 when R < 10^^ cm, we can get the corresponding 
initial isotropic energy easily as -Eiso.ini = 47r x i?ej/[27r(l — cos^^)] ~ 1.2 x 10^° ergs. We 
have already pointed out in Section 4 that our cylindrical jet looks just as powerful as an 
isotropic fireball with total kinetic energy -Bisco ~ 10^^ ergs at the GRB afterglow stage (i.e., 
when R > 10^^ cm). So our final conclusion is: if GRBs are really due to cylindrical jets, 
then an apparently energetic GRB event (eg., E^go ~ 10^^ ergs ) can be easily produced by a 
supernova- like event (i.e., -Eiso ~ 10"^*^ ergs). We even do not need to make the embarrassing 
assumption that all the kinetic energy involved in the supernova-like event is beamed into a 
narrow cone to produce the GRB, which, as noted already, is not natural. Of course, in this 
case, the GRB rate in the Universe will be much higher correspondingly. 
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Figure 1. Evolution of the Lorentz factor, 7, when v± = 0. The solid line corresponds to a cylindrical jet and the dashed line 
corresponds to a conical one. Parameters concerned here have been described in Section 4 of the main text. 





log t(s) 

Figure 3. R band afterglows from beamed GRB ejecta without lateral expansion {v± s 0). The dotted line corresponds to 
a conical jet with p = 2.5. Other lines are for cylindrical jets which differ only in the parameter p: the dashed, solid and 

dash-dotted line eorresponds to p = 3, 2.5 and 2.2 respectively. Other parameters have been given in Section 4 of the main 
text. The breaks at t ~ 10^ s in the light curves for cylindrical jets are due to cooling of electrons. 
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Figure 4. The effect of ^e, ^| and 6 on the R band hght curve (v± = case). The soHd line is plot with = 0.1, 5g = 10 ^, 
and © = 0. Each of the other lines is plot with only one parameter altered with respect to the solid line. The dashed, dash-dotted 
and dotted lines correspond to = 0.5, $g = 10"* and = 0.1 respectively. Other parameters not mentioned here have been 
described in Section 4 of the main text. 



T I I I I I I I I I I I I r 




log t(s) 

Figure 5. Evolution of the Lorentz factor, 7. The solid line corresponds to a cylindrical jet with v± = Cs- For comparison, the 
dash-dotted line corresponds to a conical jet with v± = Cg, and the dashed line is for a cylindrical jet without lateral expansion. 

Parameters are the same as in Fig. 1. 
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Figure 6. log(7 — 1) vs. logi? for jets with lateral expansion Vj_ = Cs. The solid line corresponds to a cylindrical jet and the 
dashed line corresponds to a conical one. Parameters are the same as in Fig. 1. 
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Figure 7. Schematic illustration of the lateral expansion of a cylindrical jet on the y - z plane, z-axis is the symmetry axis of 
the jet and the lateral expansion is along y-axis. Parameters are the same as those associated with the solid line in Fig. 6. 
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Figure 8. R band afterglows from beamed GRB ejecta with lateral expansion {v± = Cs). The dotted line corresponds to a 
conical jet with p = 2.5, other lines arc for cylindrical jets which differ only in the parameter p. The dash-dotted, solid and 
the dashed line corresponds to p = 2.2, 2.5 and 3.0 respectively. Other parameters not mentioned here have been described in 
Section 4 of the main text. Note that the cylindrical jets here are already non-relativistic when t > 10^ s. 
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Figure 9. The effect of ^c, and on the R band light curve {vj_ = Cs case). The solid line is plot with §o = 0.1, = 10~® 
and © = 0. Each of the other lines is plot with only one parameter altered with respect to the solid line. The dashed, dasli-dottcd 
and dotted lines correspond to = 0.5, = 10"* and = 0.1 respectively. Other parameters not mentioned here have been 
described in Section 4 of the main text. 
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